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Abstract 

The Coulomb-like backward peak of the np scattering differential 
cross section is due to one-pion exchange. Extrapolation to the pion 
pole of precise data should allow to obtain the value of the charged 
irNN coupling constant. This was classically attempted by the use of 
a smooth physical function, the Chew function, built from the cross 
section. To improve accuracy of such an extrapolation one has in- 
troduced a difference method. It consists of extrapolating the differ- 
ence between the Chew function based on experimental data and that 
built from a model where the ttNN coupling is exactly known. Here 
we cross-check to which precision can work this novel extrapolation 
method by applying it to differences between models and between 
data and models. With good reference models and for the 162 MeV 
np Uppsala single energy precise data with a normalisation error of 
2.3 %, the value of the charged nNN coupling constant is obtained 
with an accuracy close to 1.8 %. 

PACS numbers: 13.75. Cs, 13.75.Gx, 21.30.-x 

1 Introduction 



It is of crucial importance to know the precise value of the nNN coupling 
constant, both in Nuclear and in Particle physics. Together with the pion 
mass this coupling scales the nuclear interaction. Through the Goldberger- 
Treiman relation JTJ] it tests chiral symmetry. From this last relation one 
expects an accuracy in its value of about 1% as has been discussed in de- 
tails in Ref. @. Some further recent considerations on this relation and its 
implications for the nNN coupling can be found, for instance, in Ref. ||. 

In the 1980, the nNN coupling constant was believed to be well known. 
The analysis of n ± p scattering data [Q gives a value of 14.28 ± 0.18 for the 
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charged pion coupling constant. Forward dispersion relation analysis of pp 
scattering data f| led to g^o/An = 14.52 ±0.40 for the neutral pion coupling 
constant. The Nijmegen group |], 0, in the 1990's and on the basis 
of energy-dependent partial-wave analyses (PWA) of nucleon-nucleon (NN) 
scattering data, found smaller values. They obtained g^o/An = 13.47 ± 0.11 
and <7^±/47r = 13.58 ±0.05. These values were confirmed in their more 
recent NN PWA analyses |§. The Virginia Polytechnic Institute (VPI) 
group |T0| , |TTf from analysis of both ir^N and AW data has obtained also 
low values around g^/Air = 13.7. From a PWA for the 7r + p reaction a value of 
13.45(14) was recently obtained |T2]] . The very recent 7rN [IB and pion-photo- 
production Q PWA VPI analyses give values of 13.73±0.10 and 14.00±0.13 
respectively. Let us mention that some charge dependence has been also 
16|| . All the determinations which rely on the analysis of 



considered ||15l 



large data bases from a great number of experiments, with some of the data 
rejected according to certain criteria, have a very good statistical accuracy. 
It is however difficult to assert them a clear systematic uncertainty. 

Table 1: Some 7rNN coupling constants with their source. 



Reference 



Year 



Source 



Karlsruhe- Helsinki 
Kroll et al. IB 



1980 
1981 



7rp (PWA, dispersion relations) 
pp (forward dispersion relations) 



14.28(18) 
14.52(40) 



Nijmegen 



VPI 11 



Nijmegen || 
Timmermans [|l^ 



VPI 13 



VPI [HI 



1993 
1994 
1997 
1997 
1999 
1999 



pp, np (PWA) 
pp, np (PWA) 
pp, np (PWA) 

7T+ P (PWA) 

7i p (PWA, dispersion relations) 
7 p -> vrN (PWA) 



13.58(5) 

13.7 

13.54(5) 

13.45(14) 

13.73(7) 

14.00(13) 



VPI ra 



Ericson et al. 19 



1994 
1999 



GMO, vrp 
GMO, tt ± 



P 



13.75 (15) 
14.17(17) 



Uppsala 

PSI m 



Present work 



1998 
1999 
1999 



np^pn (difference method) 
np^pn (Chew/Conformal mapping 
np^pn (difference method) 



14.52(26) 
13.84(43) 
14.46(35) 



A more direct determination is the use of the Goldberger-Miyazawa- 
Oehme (GMO) sum-rule |17| which, in principle, depends directly on phys- 
ical observables. This was applied in particular in Ref. ]18| giving a value 
of g^±/47r = 13.75 ± 0.15 and very recently in Ref. |19| leading to g^±/47r = 
14.17 ± 0.17. Another direct method is based on the extrapolation to the 
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pion pole of precise data on single-energy backward differential np cross sec- 
tions. Both determinations allow a systematic discussion of statistical and 
systematic uncertainties. The use of the recent backward np Uppsala data at 
162 MeV pi, EUf and of a novel extrapolation method, the difference method, 



gives gl±/4ir = 14.52 ± 0.26 || |2T|. Fairly new analysis of the recent PSI 
backward np data, with the classical Chew extrapolation and a conformal 
mapping method leads to ^±/4vr = 13.84 ± 0.43 |22j. A summary of values 
for the coupling constant is given in Table 1. It can be seen that between the 
smallest and largest value there is a discrepancy of about 7%. The dispersion 
of the different g^±/47r values can also be judged from Fig. |TJ. 
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Kroll et al. [5] (1981) 
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VPI [18] (1994) 
Nijmegen [9] (1997) 
Timmermans [12] (1997) 
Uppsala [2] (1998) 
Pavan et al. [13] (1999) 
Arndt et al. [14] (1999) 
Ericson et al. [19] (1999). 
PSI [22] (1999) 

This work 



Figure 1: Some tcNN coupling constants as found over the last 20 years 



Let us recall that, when using the np backward data, the experimen- 
tal normalisation of the cross section is very important to the sensitivity. 
It is both the shape of the angular distribution at the most backward an- 
gles, and the absolute normalisation of the data, that are of crucial impor- 



tance 0, |20], ^T|, p23| . The Nijmegen group has strongly criticised the Uppsala 
data at 162 MeV and its extrapolated result via the difference method || |24|. 
In particular in Ref. [[| it is claimed that there is a very strong model de- 
pendence of the difference method. It is the purpose of the present study to 
check the accuracy of the difference method and in particular, when applied 
to the precise 162 MeV Uppsala data we shall demonstrate, by choosing dif- 
ferent models, that the model dependence is small and that for good reference 
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models this method leads to a precision smaller than 2%. 

The evidence that the backward peak of the np angular distribution is 
dominated by the one-pion exchange will be discussed in Sect. |2|. The de- 
termination of the kNN coupling constant through the extrapolation to the 
pion pole from models and np data is studied in Sect. [| and some conclusions 
are given in Sect. £|. 



2 Evidence for the one-pion exchange 
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Figure 2: Coulomb-like behaviour of the Svedberg Laboratory np CEX data 
2Dfl and prediction of the 'DA99' model and of other reference 



at 162 MeV 



models considered here pO . 



It was early realized that the one-pion exchange (OPE) contributes im- 
portantly to the np charge exchange (CEX) at small momentum transfer. 
In Fig. Q we have plotted the recent np CEX experimental differential cross 
section data measured at the Svedberg Laboratory at 162 MeV p|, pOf| as a 
function of q 2 , square of the momentum transfer of the neutron to the proton. 
It can be seen that there is a strong peak at very small q 2 . Note that in Fig. ^ 
q 2 is expressed in units of m,,-, the charged pion mass. This differential cross 
section has a "Coulomb like behaviour", it goes to oo, not at q 2 = as in 
the photon-exchange case, but at q 2 = — m 2 , which corresponds to the pion 
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pole of the OPE. The presence of this OPE pole very close to the physical re- 
gion led Chew in 1958 to suggest a model-independent extrapolation to this 
pole which would allow to determine the tcNN coupling constant 2"E |. 
We here refer, for the interested reader, to the very detailed and well docu- 
mented discussion concerning the OPE given in Ref . || . We shall now recall 
the different methods of extrapolation which were considered in that same 
reference. 



3 Extrapolation to the pion pole 
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Figure 3: Extrapolations to the pion pole at 162 MeV of the Chew function 
y{q 2 ) of the model DA99 |3(J for n-term polynomial fits and for the reduced 
range < q 2 < 4 m 2 of the 162 MeV np CEX Uppsala data. Here the errors 
of these experimental data have been assigned to the model DA99. 



3.1 Methods of extrapolation 

The idea to extrapolate to the pion pole is to study a smooth physical func- 
tion, the Chew function, built by multiplying the cross section by (q 2 + m 2 ) 2 . 
This removes the pole term and the extrapolation can be made more safely 
and controllably. This function can be then fitted to the data in the physical 
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region by a polynomial and extrapolated to the pole. One considers, 



. . (An) 2 sx 2 da 
y[x) = — —[x) 



d 4 



dQ 



n-1 
i=0 



Here x = q 2 + m 2 and s is the square of the total energy. At the pion pole 
x = and 

,4 



2/(0) = a = g^±/g R 
where the pseudoscalar coupling constant g 2 ±/4ir 



(2) 



14. The quantity g\ 
is a reference scale for the coupling chosen for convenience. The model- 
independent extrapolation requires accurate data with absolute normalisa- 
tion of the differential cross section. If the experimental differential cross 
section is incorrectly normalised by a factor N, the extrapolation determines 
\^Ngl±/4Tc. This is one of the most important sources of uncertainty when 
extrapolating the data. The Chew method which has been the most used in 
the past requires at least 5 terms in the q 2 expansion 0. 
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Figure 4: 5 terms Ashmore fit of the 162 MeV pseudo-data 'DA99' |3D 
the full range < q 2 < 10.1 m 2 of the np CEX Uppsala data. 



in 



A second method, which should improve the convergence of the Chew 
extrapolation, is the Ashmore method [27]. It parameterises da/dQ(x) in 
terms of a pion Born amplitude with the addition of a background term. 
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Here we use the regularised pion Born amplitudes as given in Ref. |28[]. One 
expects also an important contribution to the np CEX from the p-meson 
exchange. We then use for the Ashmore background amplitude a pole term 
with adjustable strength simulating this p-meson exchange. This expression 
is fitted to the data and gives in principle a model-independent result for the 
coupling constant. More physics is built into the procedure, so fewer terms 
should be needed. More detailed expressions can be found in Ref [fj]. 

In order to obtain an improvement in the extrapolation we have intro- 



duced the Difference Method [^TJ. It is based on the Chew function, but it 
uses the fact that an important part of the cross section behaviour is de- 
scribed by models with exactly known values for the coupling constant. It 
applies the Chew method to the difference between the function y(x) of a 
model and that of the experimental data, i.e., 



n-l 



VModelix) - y E xp(x) = ^2 d i X% - ( 3 ) 

i=0 

If gn of Eq. ([!]) is replaced by the model value g Model, one has at the pion 
pole, 

y Mod ei(o) - WO) = do = 9 ^ Mod f~ 9 - ± . (4) 

dModel 

This should decrease systematic extrapolation uncertainties and remove a 
substantial part of the non-OPE information at large momentum transfers. 
We have formally not introduced a model dependence by using such a com- 
parison function and such procedures are used in many contexts of physics to 
obtain better transparency and precision. One has to calibrate the method, 
that is, to find the precision to which the coupling constant can be deter- 
mined and the possible systematic uncertainties that are associated with the 
extrapolation procedure. 



3.2 Application to models 

We now want to cross-check to which precision can work the difference 
method in the extrapolation to the pion pole. We shall first apply it to 
models where the coupling is exactly known. This allows to investigate its 
properties and in particular its systematics. In order to determine the sys- 
tematic uncertainties in the procedures we have generated pseudo-data with 
uncertainties corresponding to the Uppsala 162 MeV experiment from 10000 
computer simulations using exact data points from different models with 
a Gaussian, random error distribution P5| . One has for a given pseudo- 



7 



Table 2: Coupling constants obtained from n-term polynomial fits to the 
31 points of several 'pseudo-data' models at 162 MeV for the reduced range 
< q 2 < 4 m 2 with the three extrapolation methods. These models are 
NN energy dependent PWA's of R. A. Arndt |30]. The unknown coupling 
constant of the model 'DA99' is to be determined while the reference models 
A13.75 and Al2.83g have gl±/4ir of 13.75 and 12.83 respectively. The column 
Sgl± is the systematic shift from the true model value. In boldface the value 
at x 2 /Ndf = 1.00, retained for the determination of gl±/4n. 



n 


X 2 /N df 




X 2 /N df 


9l±/^ 


*4± 


X 2 /N df 


^±/4vr 








'DA99' 




'Al2.83g' 






'A13.75' 












Chew Method 










4 


1.18 


11.76 ±0.38 


1.13 


10.40 ±0.43 


2.43 


1.16 


11.27 ±0.40 


2.48 


5 


1.00 


13.56 ±0.84 


1.00 


12.12 ±0.94 


0.71 


1.00 


13.04 ±0.87 


0.71 










Ashmore Method 










4 


1.01 


13.46 ±0.47 


1.01 


12.00 ±0.54 


0.83 


1.01 


13.01 ±0.50 


0.74 


5 


1.00 


13.89 ±0.94 


1.00 


12.48 ± 1.11 


0.35 


1.00 


13.40 ± 1.03 


0.35 








Difference Method 












A12.83g - 'DA99' 


A13.75 - 'DA99' 




A13.75 - 'Al2.83g 




2 


3.32 


12.15 ±0.06 


1.34 


13.51 ±0.06 




1.88 


14.15 ±0.05 


1.32 


3 


1.30 


13.09 ±0.13 


1.04 


13.85 ±0.12 




1.12 


13.62 ±0.13 


0.79 


4 


1.00 


13.97 ± 0.32 


1.00 


14.16 ± 0.32 




1.00 


13.05 ± 0.34 


0.22 



measurement m, 



with 



yFseuao-data^ = y Model^ + ^(z) (5) 



Ay m (x) = Ay u ^ sala (x)^/-2Log(R lm )cos(7TR 2m ). (6) 

In Eq. P Rim, for i=l, 2, are random numbers between and 1 when m varies 
from 1 to 10000. 

In the context of the workshop we have asked R. A. Arndt to provide us 
with different models from his NN energy dependent PWA. One of the mod- 
els, which we call A13.75, corresponds to the energy-dependent PWA of the 
pp and np data from to 400 MeV [BO] with a minimisation on gi± /An WA . 



The minimum \ 2 on the NN data is obtained for a coupling constant of 
13.75. The second model we consider is built from the previous one with all 
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parameters kept fixed, except the value of the rcNN coupling constant which 
is lowered down from 13.75 to 12.83, let us denote this model as A12.83g. 
A third model, denoted DA99, was also given to us but with an unknown 
coupling constant. We shall here apply the three methods described in sec- 
tion [3.1| in order to attempt to determine this coupling constant, considering 
the predictions of this model as pseudo-data. The differential cross section at 
162 MeV of the model DA99 is compared to that of the Uppsala experiment 
in Fig. |. 

Table 3: Same as for Table 2 but for the 54 'pseudo-data' of the full range 
< q 2 < 10.1 ml. 



n 


X 2 /N df 




X 2 /N df 


^ ± /4vr 




X 2 /N df 










'DA99' 




'Al2.83g' 






'A13.75' 












Chew Method 










5 


1.13 


12.18 ±0.29 


1.11 


10.74 ±0.33 


2.09 


1.12 


11.65 ±0.30 


2.10 


6 


1.02 


13.23 ±0.51 


1.01 


11.80 ±0.57 


1.03 


1.02 


12.71 ±0.53 


1.04 


7 


1.00 


13.86 ±0.87 


1.00 


12.42 ±0.97 


0.41 


1.00 


13.34 ±0.90 


0.41 










Ashmore Method 










4 


3.38 


11.96 ±0.33 


4.13 


9.62 ±0.41 


3.21 


3.65 


11.13 ±0.35 


2.62 


5 


1.05 


14.25 ±0.33 


1.03 


12.64 ±0.90 


0.19 


1.04 


13.69 ± 0.35 


0.06 


6 


1.03 


13.82 ±0.78 


1.02 


12.03 ±0.38 


0.80 


1.02 


12.91 ±0.84 


0.84 








Difference Method 












A12.83g - 'DA99' 


A13.75 - 'DA99' 




A13.75 - 'Al2.83g 


> 


3 


2.68 


12.23 ±0.09 


1.24 


13.55 ±0.08 




1.65 


14.14 ±0.08 


1.28 


4 


1.26 


13.31 ±0.15 


1.04 


13.93 ±0.14 




1.10 


13.48 ±0.15 


0.65 


5 


1.00 


14.07 ± 0.25 


1.00 


14.20 ± 0.25 




1.00 


12.99 ±0.27 


0.16 



Results are listed in Tables 2 and 3, n being the number of terms in the 
polynomial fit of the 'pseudo-data'. As for each calculation we performed 
10000 pseudo-experiments, \ 2 l^df is the average x 2 P er degrees of freedom, 
gl±/4:7i is the mean value of the coupling constant and the errors quoted are 
the standard deviations which, in fact, are very close to the average value of 
the error of every pseudo-experiment. We also give, for the models Al2.83g 
and A13.75, the systematic deviation Sg 2 ± of the mean value from the true 
value in the model. We have then a control on systematic extrapolation 
errors and can calibrate the corresponding corrections. The data is grouped 
in two intervals, the first one with < q 2 < 4 m\ called 'reduced range' 



9 



with 31 data points, corresponding to the previous Uppsala experiment fzl 
and the second one with < q 2 < 10.1 m 2 denoted 'full range' with 54 
data points corresponding to the latest experiment |2|, |2D[. This allows to 
examine the sensitivity and stability of the extrapolation to a given cut in 
momentum transfer and to check that it is the small q 2 region that carries an 
important part of the pion pole information. As a function of n the behaviour 
of x 2 /Ndf is characteristic: it drops quickly with increasing n to a value close 
to unity. Additional terms give only small benefits, and the data become over- 
parameterised. One can then adopt different statistical strategies leading to 
similar results. One is to take results at the minimum \ 2 ' /Njf. This minimum 
is usually a shallow one, and values of n close to n of x 2 /Ndf minimum are 
almost equally probable statistically. Another possibility is to pick up g 2 ± / 4n 
from one of the smallest values of n consistent with a \ 2 ' /Ndf well within the 
range expected from the experimental sample. We recall the reader that here 
there is about 47% probability of the experimental x 2 ' l^df to be larger than 
unity, and about 25% for it to be larger than 1.15. 

For the Chew method a good fit is performed with a fourth order polyno- 
mial in q 2 for the reduced range, but with a large systematic downward shift 
of 0.71 for both 'Al2.83g' and 'A13.75' PWA's as compared to the original 
model values. With a third order polynomial fit the statistical error becomes 
smaller, but the systematic shift is unreasonably large viz. 2.43 and 2.48 
respectively. In the full range , one or two more terms are needed to ob- 
tain a good fit. In any case a systematic shift remains even when a perfect 
fit is obtained, but at the minimum x 2 it is always less than the statistical 
and extrapolation uncertainty. The 'DA99' pseudo-data for both ranges give 
slightly different results for g 2 ±/Air at minimum x 2 , but if one applies the 
corresponding systematic shift one obtains the same value of 14. 27(86). The 
statistical and extrapolation error is rather large so we do not obtain a precise 
determination of the coupling constant using this method. Fig. ^ shows the 
fit of the model DA99 in the reduced range together with its Chew-function 
(Eq. |l]) extrapolations for n=4 (dotted line) and n=5 (solid line). The n=5 
fit is better for q 2 above 2 m 2 . 

For the Ashmore method a good description is achieved in the reduced 
range with one term less in the expansion. This shows that the physics 
beyond the 7r-exchange is reasonably described by the p-exchange, as antici- 
pated in section p.l| . The systematic shifts are similar to those of the Chew 
method. The statistical and extrapolation error, when one stands close to 
the minimum x 2 , is however smaller. Once corrected by the mean value of 
the systematic shifts of the models the pseudo-data gives a g 2 ± of 14.25(47) 
with a central value very close to the previous result. In the full range the 
needed number of terms to get a good x 2 is also smaller and the corrected 
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Difference Al2.83g - DA99 at 162 MeV 
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Figure 5: Extrapolations of the Chew function y(q 2 ) to the pion pole at 
162 MeV with the Difference Method on the model DA99 [3y| for n-term 
polynomial fits and for the reduced, a), < q 2 < 4 and full, b), < 
q 2 < 10.1 m 2 ranges. The reference function is the PWA Al2.83g f30fl . Here 
exact prediction of the model is used, the error at each point being that of 
the corresponding 162 MeV Uppsala data point 0, EU . 



value for n=5 is 14.37(35). Although the statistical and extrapolation accu- 
racy has improved, this method also appears to lack the high accuracy we 
would like to have. The excellent fit to the model DA99 is drawn in Fig. f| 
for the full range with the Ashmore parameterisation with 5 terms. 

The Difference Method should need less terms in the polynomial expan- 
sion than the two above methods, and this will give a smaller, statistical 
extrapolation error. Recall that the statistical extrapolation errors are only 
meaningful if x 2 /Ndf is close to 1. The fact that the angular distributions 
from the pseudo-data and models might be alike can help, in particular for 
large q 2 . This can add more physical information without introducing in 
principle any model dependence. We use the two reference models studied 
above. Results are given in Tables 2 and 3 for the reduced and full ranges, 
respectively. The pion-pole extrapolations of the n-term polynomial fits of 
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Figure 6: As in Fig. |5] but for the reference model 'A13.75' |30 . 



the Difference Method, on the the reduced and full ranges for the model 
DA99, are shown in Figs. |5] and |6] for the comparison models Al2.83g and 
A13.75, respectively. The error bars increase at large x, which is due to the 
multiplication of the cross section by x 2 , giving a smaller weight for the large 
q 2 region when extrapolating. The difference behaviour is slightly smoother 
with the reference model 'A13.75', however it can be seen that in all cases 4 
and 5 terms are necessary to have a good fit in the reduced and full range, 
respectively. In the full range, as one has more points, one gets a better 
statistical extrapolation error, this applies also to the difference between ref- 
erence models which gives a check on the systematic uncertainties in the 
extrapolation. Both ranges have systematic shifts below 1%. 

Averaging the values at minimum \ 2 (values in boldface in Tables 2 and 3) 
from the Difference-Method extrapolations one gets, over the reduced range, 

VNgl±/4n = 14.07 ± 0.32 (statistical + extrapolation) 
±0.11 (systematic) ± 0.16 (normalisation) 
= 14.07 ±0.37, 
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i.e. an accuracy of 2.6 %, and over the full range, 

</Ng*±/4n = 14.14 ± 0.25 (stat. + extr.) ± 0.08 (syst.) ± 0.16 (norm.) 
= 14.14 ±0.31, 

which corresponds to an accuracy of 2.2 %. The results are fairly close which 
substantiates our statement on the relevant information being nearly entirely 
at low q 2 . In view of the somewhat larger extrapolation uncertainty in the 
case of the reduced range, we take the full range value, 14.14(31), which 
compares quite well with the exact value of the model DA99 which is 14.28. 
Using the experimentally given statistics of the 162 MeV Uppsala data, the 
difference method has allowed us to determine the unknown DA99 coupling 
to less than 1 % within an uncertainty of 2.2 %. 

3.3 Application to data 



Table 4: g 2 ±/4ir and f 2 /^ of the PWA comparison models |30| used here 
together with their total y 2 and x 2 /data on the 3747 np data below 400 
MeV. 



Model 



A12.83g A12.83 A13.75 A14.28 Al4.28g 



xl P (total) 
X 2 np / data 



12.83 


12.83 


13.75 


14.28 


14.28 


0.071 


0.071 


0.076 


0.079 


0.079 


9282 


5087 


4975 


4958 


6149 


2.48 


1.36 


1.33 


1.32 


1.64 



We shall now apply the difference method to the 162 MeV Uppsala 
data g H to cross-check the precision of the g^±/47r determination of 



Ref. 0. We here consider different reference models than those used in 0. 
Besides the three comparison models considered in the previous section we 
also use two more PWA fit of the pp and np data from to 400 MeV with 
^±/4tt fixed to 12.83 (model A12.83) and 14.28 (model A14.28) 0. In 
Table 4 we give the total x 2 an d the x 2 /data on the 3747 np data consid- 
ered for all the models we study here. We also recall their g 2 ±/47r together 
with the corresponding pseudovector coupling f 2 /4:7r, related to g^^/Air by 
g 2 ± = f 2 (2M MN /m n ±) 2 with M N average proton- neutron mass. As for the 
model Al2.83g, the model Al4.28g (DA99) was obtained from the model 
A13.75 with all parameters kept fixed but increasing the 7rNN coupling from 
13.75 to 14.28. Predictions of all these comparison models for the 162 MeV 
Uppsala data are shown in Fig. 0. 
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Figure 7: As in Fig. |5| but for the difference between model Al2.83g |30| and 
Uppsala data fl2L 120]. 



Results for the reduced and full range are listed in Tables 5 and 6 respec- 
tively. The A12.83g comparison model requires 4 terms in the reduced range 
and 5 in the full one. Figure [7] shows that there is some edge effect in the re- 
duced range at q 2 = 4m 2 . For the model A12.83 one needs also 4 terms in the 
reduced range. The obtained value, 14.37(31) is quite consistent with those 
determined in the model Al2.83g, either in the reduced range, 14.69(31) or 
in the full range, 14.41(24). In the full range the minimum is very shallow 
and the n = 4 (x 2 /N df = 1.18) and n = 5 (x 2 /N df = 1.16) values, 13.65(14) 
and 13.94(25), respectively, are compatible. The lower value, 13.65(14), is 
however not compatible with the value, 14.37(31) of the reduced range at 
X 2 /Ndf minimum, We shall then retain, in the full range for the A12.83 
model, the n = 5 determination, 13.94(25). Figure |8] shows as before the 
necessity to have at least 4 or 5 terms to obtain a good fit. In the reduced 
range 3 terms are here sufficient with the reference model A13.75 leading to 
g 2 ±/4ir = 14.50(12) while in the full range x 2 minimum is reached with 4 
terms with a value of 14.38(14). Figure |9] shows in both range a smoother 
behaviour than before. In the A14.28 and Al4.28g case one needs the same 
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Table 5: As in Table 2, but, i) for the difference method only, ii) the 31 points 
of the pseudo-data of the model DA99 are replaced by those of the 162 MeV 
Uppsala data 0, |2tJ, iii) with 3 more Arndt PWA reference models [ |30[] , viz., 



A12.83, A14.28 and Al4.28g (same as 'DA99'). Models A12.83 and A14.28 
(as also A14.28g) have gl±/4ir of 12.83 and 14.28, respectively. In boldface 
the experimental values retained for the determination of g^±/A.ir. 









Difference Method 








n 


X 2 /N df 




X 2 /N df 






X 2 /N df 








A12.83j 


I — Uppsala 


A12.83 


— Uppsala 




A13.75 


— Uppsala 




2 


1.60 


13.45 ±0.05 


1.13 


13.85 ±0.05 




1.13 


14.69 ±0.05 




3 


1.34 


13.79 ±0.13 


1.15 


13.75 ±0.13 




1.06 


14.50 ± 0.12 




4 


1.03 


14.69 ± 0.31 


1.03 


14.37 ± 0.31 




1.03 


14.87 ±0.30 




5 


1.07 


14.50 ±0.79 


1.07 


14.30 ±0.80 




1.06 


14.54 ± 0.78 






A14.28 


— Uppsala 


A14.28* 


? — Uppsala 




A13.75 - '12.83' 




2 


1.17 


15.20 ±0.05 


1.72 


15.42 ±0.05 




1.05 


13.76 ±0.05 


0.93 


3 


1.02 


14.96 ± 0.12 


1.00 


14.94 ± 0.12 




1.03 


13.65 ±0.13 


0.82 


4 


1.04 


15.19 ±0.30 


1.03 


15.00 ± 0.30 




1.01 


13.40 ± 0.33 


0.57 


5 


1.06 


14.72 ± 0.77 


1.06 


14.59 ± 0.78 




1.00 


13.09 ±0.87 


0.26 








A12.83 - 'A14.28' 




A13.75 - 'A14.28' 




2 






1.14 


12.86 ±0.06 


-1.42 


1.02 


13.76 ±0.05 


-0.52 


3 






1.07 


13.04 ±0.13 


-1.24 


1.01 


13.83 ±0.12 


-0.45 


4 






1.01 


13.43 ± 0.33 


-0.85 


1.00 


13.96 ±0.32 


-0.32 


5 






1.00 


13.83 ± 0.83 


-0.45 


1.00 


14.09 ±0.81 


-0.19 



number of terms as for A13.75 but the corresponding gr^ ± /47r are somewhat 
larger and not always compatible with previous values. Curves for the ex- 
trapolation are shown in Figs. [H] and [11]. Applying the difference method 
between models (see also Tables 2 and 3) shows relatively large systematics 
which allows to understand this dispersion. 

Let us here remind the reader that the difference method analysis is quite 
consistent. One can check, in the Tables 2, 3, 5 and 6, that, at a given n 
value, the difference between the results of the difference method for the 
reference models and data is very close to the systematic shift between the 
two models. One has, 

s£±/47r(ra, Model A - Uppsala) - g£±/47r(n, Model B - Uppsala) 

~ 5gl±(n, Model A - Model B). (7) 

For instance if we compare in the full range the n = 4, (A13.75 — Uppsala) 
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Table 6: As in Table 5 but for the 54 Uppsala data points at 162 Mev of the 
full range < q 2 < 10.1 m 2 . 
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11 


,.2 / AT 


9 n ±/ 47T 


X 2 /N df 






X 2 /N df 




An 2 




AlZ.oof 


5 — Uppsala 


A12.83 


— Uppsala 




A13.75 


— Uppsala 




3 


1.30 


13.65 ± 0.08 


1.36 


14.02 ±0.08 




1.43 


14.84 ±0.07 




4 


1.30 


13.79 ±0.14 


1.18 


13.65 ±0.14 




1.12 


14.38 ± 0.14 




5 


1.14 


14.41 ± 0.24 


1.16 


13.94 ± 0.25 




1.14 


14.53 ± 0.24 




6 


1.15 


14.74 ±0.46 


1.16 


14.42 ± 0.47 




1.15 


14.80 ± 0.45 






A14.28 


— Uppsala 


A14.28j 


I ~ Uppsala 




A13.75 - '12.83' 




3 


1.51 


15.33 ±0.07 


2.15 


15.53 ±0.07 




1.03 


13.74 ±0.08 


0.91 


4 


1.11 


14.83 ± 0.13 


1.12 


14.75 ± 0.13 




1.02 


13.63 ±0.14 


0.80 


5 


1.13 


14.89 ± 0.24 


1.13 


14.62 ± 0.24 




1.01 


13.49 ± 0.26 


0.66 


6 


1.15 


15.05 ±0.45 


1.15 


14.86 ± 0.45 




1.00 


13.27 ±0.50 


0.44 








A12.83 - 'A14.28' 




A13.75 - 'A14.28' 




3 






1.08 


12.89 ±0.09 


-1.41 


1.01 


13.77 ±0.08 


-0.51 


4 






1.04 


13.07 ±0.15 


-1.21 


1.01 


13.83 ±0.14 


-0.45 


5 






1.02 


13.30 ±0.26 


-1.13 


1.00 


13.92 ±0.25 


-0.36 


6 






1.01 


13.62 ± 0.49 


-0.66 


1.00 


14.02 ± 0.48 


-0.26 



result to that of (Al2.83g — Uppsala) we have a difference (see Table 6) of 
0.59 (14.38 - 13.79) which is close to the systematic shift between these models 
.65 (see Table 3 ). In Ref. [0] where the comparison models were the Nijmegen 
potential |31j, the Nijmegen f32[ (NI93) and Virginia |3|, (SM95) energy 



dependent PWA's, dispersion of results were smaller. This could be traced to 
the fact that these models, where g 2 ±/47r has been minimised with respect to 
the A^ A" data, have a high q 2 momentum more similar to that of the Uppsala 
data. 

Summarising these results we take for g 2 the following average in the 
reduced range 

gl ± /4n = -[14.69(31) + 14.37(31) + 14.50(12) + 14.96(12) + 14.94(12)] 

o 

= 14.69(20) (8) 

It is here only the second value (from A12. 83— Uppsala) which is slightly 
outside the range of the 2 last values from A14.28 and A14.28g 

Taking half of the average of |5<7^±| between models at x 2 /^df — 1-00 for 
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Figure 8: As in Fig. ||] but for the difference between model A12.83 and 



Uppsala data pi 20 . 



the estimation of the systematic uncertainty (see Table 5) we obtain 

gl ± /A7r = 14.69 ± 0.20 (stat. + extr.) ± 0.15 (syst.) ± 0.17 (norm.) 

= 14.69(30) (9) 

i.e. an accuracy of 2 %. For the full range (see Table 6) we have 

gl±/4x = -[14.41(24) + 13.94(25) + 14.38(14) + 14.83(13) + 14.75(13)1 
5 



14.46(18) 



(10) 



Adding the systematic (estimated as above) and normalisation errors we 
obtain 



gl ± /4n = 14.46 ± 0.18 (stat.+ extr.) ± 0.15 (syst.) ± 0.17 (norm.) 
= 14.46(29) 



IF 
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Figure 9: As in Fig. |5] but for the difference between model A13.75 |K| and 



Uppsala data || |20 . 



i.e. again an accuracy of 2%. Note that for A12.83 — A14.28 we need to go 
to n = 7 to get x 2 /N df = 1.00 with a 6g*± of -0.36. The value of Eq. QTTJ) 
is to be compared to the value we determined in Ref. @, viz. 



14.52 ± 0.13 (stat.+ extr.) ± 0.15 (syst.) ± 0.17 (norm.) 
14.52(26) (12) 



i.e. an accuracy of 1.8%. It is seen that both determination, Eqs. ( i~l~D 
and (|12"D, are very close. The determination of Ref. 0, Eq. ( O ) has a better 
statistical extrapolation error which can be understood as the comparison 
model have possibly a better determined high q 2 behaviour as just mentioned 
above. 



4 Conclusions 

This analysis of the 162 Mev Uppsala precise experiment on np charge ex- 
change demonstrates here again that such data can be used for a direct and 
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Figure 10: As in Fig. |5] but for the difference between model A14.28 |30| and 



Uppsala data M, 20 . 



accurate determination of the irNN coupling constant. We reproduce the 
original coupling constant using the present procedures for pseudo-data from 
PWA's as exemplified with the model DA99. The extrapolation error in- 
creases with the number of parameters. Our value is 7% larger than the 
Nijmegen || result g\± = 13. 58 ±0.05, but it is consistent with values given 
in earlier data compilations based on the analysis of ttN and AW scattering 



data M5 



The data have been used to determine a precise value for the charged 
ttNN coupling constant using extrapolation to the pion pole. Using the most 
accurate extrapolation method, the Difference Method, we find \fN gt± = 
14.46 ± 0.18 (ff± = 0.0800 ± 0.001) with a systematic error of about ±0.15 
(±0.0008) and a normalisation uncertainty of ±0.17 (±0.0009). We do repro- 
duce the input coupling constants of models using equivalent pseudo-data. 
The practical usefulness of the method, its precision and its relative insensi- 
tivity to systematics appear to be under control. The pseudo-data demon- 
strate that considerable precision is achieved statistically at a single energy. 
The absolute normalisation of the data is nevertheless crucial. The precision 
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Figure 11: As in Fig. |5| but for the difference between model Al4.28g [30 
and Uppsala data P, p0| . 



of the method used here has not yet reached its theoretical limit, but we can 
point out the key information necessary for this in the AW sector. We need 
as precise as possible unpolarised differential cross sections with an absolute 
normalisation of 1 to 2% to reach a precision of about 1% in the coupling 
constant. For the Uppsala data an accurate normalisation 0, |(J of 2.3 % 
was obtained using integration over the angular distribution. It was per- 
formed on the part they measured (from 72° to 180°) and on the remaining 
one calculated from PWA's and models. The result was scaled to the well 
known experimentally total cross section [RBI. It is important to extend the 



angular range of data, to be able to achieve an improved normalisation. 

For that purpose there is, as we have heard [OT in this workshop, an 
experiment in progress at The Svedberg Laboratory. It will measure the np 
differential cross section in the forward hemisphere. This should hopefully 
determine the normalisation to 1 % and then lower down the accuracy on 
<7^±/47r from 1.8 % to 1.5 %. To the question which was asked at this work- 
shop "Does the difference method work at the 1 % level?" the answer is so 
far no, but as we have demonstrate it can work to less than 2 %, in particu- 
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lar with good reference models we obtained a precision of 1.8 %. Here with 
different reference models the precision reached was 2 %. Contrary to what 
was claimed in Ref. M we do not see a large model dependence, even using 
"extreme" models as the Al2.83g with ^±/47T=12.83 or / c 2 /4tt=0.071. The 
model dependence is relatively small as can be judged from a comparison 
of the result here and that of Ref. where while using different reference 
models the results agree within less than 0.5 %. 

Let us also mention, as we were told [|37| in this workshop that np 
backward measurement with tagged neutron beams (which leads to abso- 
lute normalisation) are in progress at IUCF for 185 > < 195 MeV and 
90° < 9 < 180°. If the very backward steeper shape of the Uppsala data as 
compared with earlier data [fj, [21], is confirmed together with its absolute 
normalisation then the itNN coupling constant cannot be as low as found, for 
instance, by the Nijmegen group. If the Uppsala data is correct and if the 
coupling constant is small then to reconciliate both, either the total cross 
section could be off or the angular distribution in the forward hemisphere 
could be different of what it is believed so far from PWA analyses. Hopefully 
the two above mentioned experiments in progress should help to clarify the 
situation. There will be also pp and np spin-transfer measurements which 



could be used to precise the itNN coupling constant [38 . 

In principle, an experiment at one single energy is enough to determine 
<j^±/47r, for all energies contain similar information. Although the method 
of analysis seems to work well, it is, however, useful to deduce the coupling 
constant from data at several energies. This would increase confidence that 
some unexpected systematic effect influences the conclusion. 
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